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ABSTRACT 


The  dynamic  equations  for  a  Liquid  Metal  Fast  Breeder 
Nuclear  Reactor  (LMFBR)  Power  Plant  are  derived.   These  non- 
linear differential  equations  with  typical  design  values  for 
the  system  parameters  are  solved  on  a  COMCOR  Ci  5000  analog 
computer  for  the  characteristic  transients  of  the  basic  plant 
and  the  results  compared  with  a  recently  proposed  LMFBR  design 
A  digital  computer  automatic  control  program  is  developed 
which  uses,  in  addition  to  the  ideas  of  derivative  and  pro- 
portional control,  precomputed  control  trajectories.   The 
digital  controller  design  is  implemented  by  applying  FORTRAN 
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tiveness  of  the  strategy  for  the  automatic  controller  is 
demonstrated  by  using  various  step  changes  in  the  load 
demand. 

A  special  software  routine  was  written  to  automatically 
scale  the  analog  computer  and  proved  to  be  quite  useful  in 
improving  the  operating  efficiency  of  the  hybrid  computer 
simulation. 
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I.   INTRODUCTION 

In  view  of  the  enormous  increase  of  world  power  consumption 
and  the  fast  decrease  of  the  resources  which  may  be  used  to 
satisfy  this  growing  power  demand,  it  became  necessary  to 
develop  and  build  nuclear  power  plants.   The  reactors  used  to 
generated  the  power  for  these  plants  are  up  to  now  almost 
exclusively  thermal  reactors  using  uranium  as  fissionable 
material.   There  are  3  fissionable  isotopes:   U  (235),  Pu  (239), 
and  U  (2  33) .   U  (2  35)  is  found  in  natural  uranium  and  can  be 
used  in  this  state  or  in  an  enriched  form,  i.e.,  its  percentage 
is  increased  beyond  that  in  natural  uranium.   Pu  (239)  is 
manufactured  by  nuclear  reaction  from  U  (238)  .   U  (233)  is 
similarily  manufactored  from  Th  (2  32)  ,  an  isotope  which  is  also 
found  in  natural  uranium. 

The  worlds  supply  of  nuclear  fuel  consists  of  25x10   tons 
of  natural  uranium,  mostly  U  (2  38) ,  and  1x10   tons  of  thorium, 
mostly  Th  (232)  [10].   The  natural  uranium  consists  of  99.29% 
U  (238)  and  0.71%  U  (235).   Since  U  (235)  is  the  only  naturally 
fissionable  isotope  the  use  of  nuclear  fuels  would  be  limited 
to  0.71%  of  25x10   tons.   This  means  that,  using  only  the 
present  thermal  type  reactors  with  their  inefficient  nuclear 
economy  a  shortage  of  nuclear  fuel  is  lightly  to  occur  by  the 
year  2000. 

An  increased  nuclear  economy  in  an  uranium  reactor  can  be 
achieved  by  using  fuel  breeding.   In  this  process  the  nonfiss: 


capture  of  neutrons  by  a  fertile  U  C238)  nucleus  results  in 
the  generation  of  fissionable  Pu  (239) .   Fertile  Th  (232)  can 
also  be  converted  to  fissionable  U  (2  33) .   If  more  fissionable 
nuclei  are  produced  than  are  consumed  the  reactor  in  which  this 
process  takes  place  is  called  a  Breeder  Reactor.   If  the  fission 
process  depends  on  fast  neutrons  as  in  the  above  described 
processes  the  reactor  is  a  Fast  Breeder  Reactor  (FBR) . 

The  control  Problem  of  a  FBR  is  particularly  acute  because 
of  two  important  facts.   First  of  all  due  to  the  high  energies 
required  of  the  neutrons  in  order  to  produce  fission,  the  power 
per  unit  volume  is  much  larger  than  in  the  presently  used 
thermal  type  reactors  [10] .   The  high  power  density  results  in 
steep  temperature  gradients,  large  coolant  temperature  rises 
and  therefore  coolant  distribution  problems. 

The  second  important  fact  is  the  lack  of  inherent  load 
following  capability  of  this  type  of  plant  [X]  .   This  is  due 
to  the  fact  that  the  load  fluctuations  do  not  feedback  as 
efficiently  through  the  coolant  medium  as  they  do  in  thermal 
type  reactors .   For  these  reasons  constant  automatic  super- 
vision and  control  of  a  fast  breeder  reactor  is  necessary. 


To  design  a  controller  for  a  nuclear  power  plant  it  is 
necessary  first  to  obtain  a  mathematical  model.   In  this  report, 
the  dynamic  equation  of  a  LMFBR  are  derived.   The  next  step  in 
the  design  procedure  was  to  determine  the  dynamic  behavior 
of  the  complete  plant  following  selected  changes  in  both  the 
controlling  inputs  (reactivity,  valve  opening)  and  in  the  load. 
Because  of  the  nonlinearity  of  the  system,  model  simulation 
was  required  in  order  to  obtain  a  solution  for  the  dynamic 


response.  From  this  simulation  conclusions  were  drawn  about 
the  characteristics  which  the  controller  must  have  and  about 
the  required  control  strategy. 

The  real-time  solution  of  the  dynamic  equations  is  performed 
by  the  analog  computer  CCOMCOR  Ci  5000)  in  a  hybrid  linkage 
with  the  digital  computer  CXDS  9  300) . 

The  controller  design  developed  in  this  report  uses  the 
digital  computer  as  the  controlling  element.   The  economic 
advantage  of  using  a  digital  controller  results  from  the  pros- 
pect of  operating  a  general  purpose  type  computer  in  the  time- 
share  mode  to  perform  both  the  control  decisions  as  well  as  the 
monitoring  data  processing  and  such  other  related  functions 
that  are  presently  being  performed  by  such  digital  computers. 

A  digital  computer  automatic  control  program  is  then 
developed  on  the  basis  of  the  results  of  the  simulated  plauL 
response.   The  design  strategy  includes  in  addition  to  the 
ideas  of  derivative  and  proportional  control ,  precomputed 
control  trajectories.   The  digital  controller  design  was 
implemented  by  applying  FORTRAN  IV  programming  to  the  XDS  9  300 
digital  computer.   The  effectiveness  of  the  strategy  for  the 
automatic  controller  is  demonstrated  by  using  various  step 
changes  in  the  load  demand. 

A  special  software  routine  was  written  to  automatically 
scale  the  analog  computer  and  proved  to  be  quite  useful  in 
improving  the  operating  efficiency  of  the  hybrid  computer 
simulation. 
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II.-  PLANT  DESCRIPTION 


The  plant  considered  in  this  investigation  consists  of: 

1)  Fast  Breeder  Nuclear  reactor  with  control  rods 

2)  Primary  sodium  coolant  loop 

3)  Heat  exchanger  (counterf low) 

4)  Secondary  sodium  coolant  loop  bypass  and  bypass 
control  valve 

5)  Steam  generator 

6)  Steam  loop  with  turbine  and  condensor 
Additional  components  which  are  found  in  nuclear  reactor 

.power  plants  such  as  pumps,  intermediate  heat  exchangers, 


jjiGssuic  uanjcs  ,  ecj uaJ- j- z er  t_anj».3  ,  s  <-orage  <_an..v.3  , 


superheater,  and  intermediate  turbine  stages,  were  not  included 
since  they  would  not.  have  served  the  purpose  of  this  investi- 
gation.  By  adding  these  parts,  the  basic  structure  of  the 
system  would  not  be  changed  and  would  still  be  valid. 
A  flow  diagram  of  the  plant  is  given  below. 
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R  =  Fast  Breeder  nuclear  reactor 

CR  =  Control  rods 

HX  =  Heat  exchanger 

SGen  =  Steam  generator 

T  =  Turbine 

C  =  Condenser 

Nal  =  Primary  sodium  loop 

Na2  =  Secondary  sodium  loop 

HpO  =  Steam  circuit 

B  =  Sodium  bypass  valve 

TF  -  Average  fuel  temperature 

TAV  =  Average  coolant  temperature  in  reactor  core 

TCI  =  Reactor  inlet  temperature 

TH1  =  Reactor  outlet  temperature 

TBIl  =  Inlet   temperature  heat  exchanger  primary  side 

TB01  =  Outlet  temperature  heat  exchanger  primary  side 

TBI  =  Average  temperature  heat  exchanger  primary  side 

TH2  -  Outlet  temperature  heat  exchanger  secondary  side 

TC2  -  Inlet  temperature  heat  exchanger  secondary  side 

TB2  =  Average  temperature  heat  exchanger  secondary  side 

TBI2  =  Inlet  temperature  steam  generator  primary  side 

TB02  =  Outlet  temperature  steam  generator  primary  side 

TB3  =  Average  temperature  steam  generator  primary  side 

TS  =  Average  temperature  steam  generator  secondary 
side 
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III.   MATHEMATICAL  MODEL  -  PLANT  DYNAMICS 


A.   REACTOR  CORE 

1.   Neutron  Kinetics 


The  neutron  kinetic  equation  with  one  group  approxima- 
tion [1]  is  given  as 


dnO)_  =  3KCt)-3  n(t)  +  x  c(t)  (m>1) 


^1^1  =  B  n(t)  _  X  c(t)  (III. 2) 

at      x, 

where    n(t)  =  neutron  density 
6K(t)  =  reactivity 

3  =  fraction  of  delayed  neutrons 
I   =   neutron  lifetime 
X   =  average  decay  constant 
c(t)  =  average  concentration  of  delayed  neutrons 
The  steady  state  conditions  are  given  by 

6k  (t)  =  0 


^n   =  Ac 
e   o     o 


where     n   =  steady  state  neutron  density 

c   =  average  steady  state  concentration  of 
delayed  neutrons 
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When  setting 


•^i^-  =  K(t)  . 

p 

=  reactivity  in  dollar  units 
n(t) 


n 
o 


-   NCt) 


cCti=    D(t) 


c 
o 


the    following   equations    are   obtained 


Sj|£L  =    |    CK(t)     N(t)     -    N(t)     +    D(t)  (III. 3) 


d^t}    ■    A     (N(t)     -    D(t))  (III. 4) 

The  numerical  values  for  a  fast  reactor  for  6  and  £  are  approx- 
imately given  as 

B  s  3  •  10~3 
I    =    6  •  10"7 

4 


Therefore  j   -  .5  •  10 

If  the  differential  equation  for  N(t)  is  brought  into  the  form 

I    dN(t) 


0  =  K(t)  N(t)  -  N(t)  +  D(t) 


3   dt 


The  last  term,  being  small  compared  to  the  other  three  terms 
may  be  neglected.  Using  this  approximation  the  equation  for 
the  neutron  kinetics  in  the  reactor  core  are  then  given  as 

^~^-  =  X    (N(t)  -  D(t))  (III. 5) 

dt 

N;t)  -  mf  (iii-6) 
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with  the  steady  state,  initial  conditions 

N(0)  =  D(0)  =  1  (III. 7) 

2 .   Reactivity  Coefficients 

In  addition  to  the  controlled  reactivity  input  applied 
by  positioning  of  the  control  rods,  a  nuclear  reactor  has 
inherent  reactivity  feedback,  produced  by  thermohydraulic 
effects.   The  total  reactivity  is  given  as 

K  (t)  =  KR(t)  +  Kp(t)  (III. 8) 

where    K  (t)  =  controlled  reactivity  input 

KD(t)  =  inherent  reactivity  feedback  coefficients 

r 

In  a  fast  reactor  the  following  reactivity  coefficients  are 
of  importance 

a)  Doppler  coefficient 

b)  Fuel  axial  expansion  coefficient 

c)  Fuel  radial  expansion  coefficient 

d)  Clad  temperature  coefficient 

e)  Sodium  temperature  coefficient 

a  block  diagram  for  the  stability  model  is  given  in  Fig.  2. 
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Fig.  2.  Block  Diagram  for  reactivity  feed  back  paths 


The  numerical  values  of  the  reactivity  coefficients  are  given 
as  [2] 

a)  Doppler  coefficient        -  1.39    cents/  %  power 

b)  Fuel  axial  expansions       -   .32    cents/  %  power 
coefficient 

c)  Clad  temperature  +   .0138  cents/  %  power 
coefficient 
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d)  Sodium  temperature         +   .0057  cents/  %  power 
coefficient 

e)  Radial  expansion  -   .0043  cents/  %  power 
coefficient 

From  the  above  data  it  is  seen  that  coefficients  c,  d,  and  e 
are  several  orders  of  magnitude  smaller  than  a  and  b.   There- 
fore, these  reactivity  coefficients  will  be  neglected  and  it 
will  be  shown  in  section  IV  that  this  does  not  lead  to  any 
noticeable  inaccuracy.   The  reactivity  feedback  coefficient 
may  now  be  expressed  as 

KF  "  KD  +  KEXP  '  '     (III-9) 

where    K   =  doppler  coefficient 

K„.T_  =  fuel  axial  exoansion  coetticient 
and  the  equation  for  the  total  reactivity  is  given  as 

KT  =KR  1-  KD  +  KExp  (III. 10) 

If  the  numerical  values  for  the  reactivity  coefficients  are 
introduced  along  with  the  dollar  units  of  reactivity  (see 
Appendix  A) ,  the  final  equation  for  the  reactivity  feedback 
is  given  as 

Km  =  KD  -  .0014  TF  (III. 11) 

with  the  control  block  diagram  given  in  Fig.  3. 
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Fig.  3.   Reactivity  Feedback  Model 
(Simplified) . 


3 .   Heat  Generation 

The  rate  of  heat  produced  per  unit  volume  of  fuel  is 
called  the  volumetric  thermal  source  strength  0,,  and  is  given 
by  [4] 


-N 


n(t) 


* 


(III. 12) 


where    G   =  energy  per  fission,  Mev 

3 
n(t)   =  number  of  fissionable  fuel  nuclei/cm 

a    =  Microscopic  fission  cross  section  of  the  fuel 

2 

$   =  neutron  flux  per  cm   sec. 

If  n  is  taken  as  a  constant  for  a  certain  fuel  element,  it  is 
seen  from  the  above  relation,  that  the  generated  heat  is 
proportional  to  the  product  of  the  fission  cross  sec 
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and  the  neutron  flux  cj>.   The  neutron  flux  Is  dependent  upon 
the  fissionable  species  used  and  the  energy  of  fission  neutrons, 
which  in  turn,  is  dependent  upon  the  reactor  temperature.   If 
o^   if  determined  for  a  certain  reactor,  it  is  seen  that  the 
heat  generated  at  any  point  is  directly  proportional  to  the 
neutron  flux  at  that  point,  which  is  directly  proportional  to 
the  earlier  used  N  (t)  .   The  generated  heat  may  therefore  be 
written  as 

QN  =  A  N(t)  (III. 13) 

where    A  =  proportionality  factor,  which  is  constant  if 
Or;   is  assumed  to  be  a  constant 

B.   HEAT  TRANSFER 
1.   Reactor 

a.   Fuel  to  Coolant 

The  core  of  a  fast  nuclear  reactor  consists  of  a 
great  number  of  fuel  elements.   In  a  proposed  1  MW  Fast 
Breeder  Reactor  [2]  for  example,  there  are  4  70  fueld  rods,  each 
fuel  section  (core  length)  is  2  feet  long  and  the  rod  has  an 
overall  diameter  of  .25  inches.   The  fuel  is  enclosed  in  a 
cladding  of  wall  thickness  of  .015  inches. 

Due   to  the  size  of  the  fuel  element  the  radial 
flux  distribution  within  the  core  may  be  neglected.   The  axial 
flux  distribution  is  shown  in  Figure  4. 
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Fig.  4.   Reactore  Core  and  Axial 
Flux  Distribution. 


If  the  axial  variation  in  flux  is  taken  into  account,  it  is 
seen  from  figure  5  that  the  generated  heat  and  therefore  the 
temperature  in  the  fuel  element,  cladding  and  the  coolant 


varies  in  the  axial  direction 
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Fuel 
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Fig.  5.   Temperature  variation  along  fuel 
element  (center,  cladding, 
coolant) . 
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In  this  investigation,  however,  since  the  reactor 
is  only  a  part  of  the  total  power  plant,  these  special  effects 
will  be  neglected.   Only  average  temperatures  will  be 
considered. 


Fig.  6.   Temperature  profile  in 
fuel  element. 


Figure  6  shows  a  typical  temperature  profile  within 
the  fuel  element,  cladding  and  coolant  together  with  the  average 
temperature  profile  considered  here.   The  effect  of  the  clad- 
ding is  neglected,  which  is  justified  in  view  of  the  earlier 
given  numerical  values.   Taking  the  average  temperatures 
simplifies  the  mathematical  solution  without  introducing 
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considerable  errors.   If  more  accuracy  is  desired,  the  obtained 
temperature  values  could  always  be  corrected  with  a  multipli- 
cation factor  for  point  values  or  for  optimal  accuracy  the 
temperature  distribution  model  may  be  obtained  by  the  solution 
of  partial  differential  equations  [3]. 

Under  the  above  assumptions  the  heat  balance  for 
a  single  fuel  element,  whose  diameter  is  small  with  respect  to 
the  core,  may  be  written  for  the  steady  state  heat  transfer 
as  : 


Q„  2-nr   dr  L  =  Q  ,  ,   -  Q 
*F  *r+dr    r 


(III. 14) 


Fig.  7.    Cross  section  of  cylindrial 
fuel  element. 
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The  solution  to  this  equation  with  the  boundary  condition 


dT 

dr 


r=0 


r=0 


=  TF 


is  (see  Appendix  B) 


QF  =  K  A  (TF  -  TAV) 


(III. 15) 


where     K  =  Thermal  conductivity,  in  this  case  an  average 
value  for  the  fuel  cladding  and  coolant 

cal 


sec  cm   °C 

A  =  surface  area  of  the  fuel  element  [cmz] 
The  thermal  conductivity  of  the  fuel  and  cladding  and  the 
physical  properties  of  the  coolant  (density,  viscosity,  and 
specific  heat)  are  considered  to  be  constant  and  independent 
of  temperature,  which  is  a  valid  assumption  for  medium  power 
variations  of  the  reactor.   For  the  above  assumptions  the 
product  K  .  A  may  be  considered  constant,  depending  only  upon 
the  physical  size  of  the  fuel  element.   For  the  dynamic  case 
of  the  heat  transfer  the  heat  capacitance  of  the  fuel  must  be 
taken  into  account.   Using  the  analogy  of  an  electric  circuit 
for  heat  flow,  where  heat  flow  is  equivalent  to  current, 
temperature  equivalent  to  voltage,  and  heat  sink  (dissipation) 
equivalent  to  a  resistor,  the  following  relationship  can  be 
written 


dTF 

N     f  ,.     *F 

dt 


(III. IS) 
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Fig.  8.  Heat  flow  model. 


The  definition  of  the  heat  capacitance  is  [7] 


C  = 


A  T 


cal 
sec 


Therefore  the  above  equation  may  be  written  as 


~-  =   ~  (TAV  -  TF)  +  A'  N 


(III. 17) 


u                  1    K  A 
where    —  =  


cal 


2 

sec  cm   °C 


cm 


sec 


cal 


=  fuel  time  constant 
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A' 


A_ 
c, 


=  Proportionality  factor 
b.   Heat  Removal 

The  amount  of  heat  removed  by  the  coolant  is  given 


as 


Q  =  3  •  A, 


V 


AT 


(III. 18) 


where    3  =  density  of  coolant 


gm 


cm 


A 


cross  sectional  area  of  coolant  passage  [cm  ] 


V  =  coolant  speed 


cm 
sec 


C  =  specific  heat  of  coolant 


cal 


gm  °C 


J 


AT  =  coolant  temperature  rise  within  the  element  [°C] 
As  shown  earlier  the  temperature  rise  of  the  coolant  is  not 
a  linear  function  of  the  position  of  a  unit  element  within 
the  coolant  channel.   A  second  nonlinearity  accours  in  the 
coolant  channel  which  is  shown  in  Fig.  9.   The  temperature  of 


Fuel 


Fuel 
Fig.  9.    Temperature  in  coolant  path 
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coolant  at  the  walls  of  the  fuel  element  are  higher  than  at  the 
center  of  the  coolant  channel.   This  nonlinearity  may  be 
neglected  for  the  reasons  stated  earlier,  and  an  average  tem- 
perature profile  will  be  assumed  as  shown  in  Fig.  9.   The 
specific  heat  of  the  coolant  and  the  coolant  density  will  be 
assumed  constant.   The  heat  balance  equation  for  the  dynamic 
case  is  again  derived  with  the  aid  of  an  electrical  equivalent 
as  shown  in  Fig.  10. 
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Fig.  10.  Electrical  equivalent  for  heat  removal 
by  coolant  from  core . 


The  heat  flow  equation  is 


q   ,  c   *™  +  Q 
WF     c   dt     yAV 


1.19) 
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where   QAV  =  Q    (.from  equation  m.  18) 


which  may  also  be  written  as 


Q 


where 


AV 


m  = 


m  c  AT 

m  c  (TCI  -  TH1) 

coolant  flow  rate 


gm 


sec 


After  rearranging  and  using  the  relation 


TAV  = 


TCI  +  TH1 


to  eliminate  TH1,  the  following  equation  is  obtained 


dTAV 
dt 


+  —V   TAV  +  —  TF  +  ^-r  TCI 

T2    T2  T2       T2 


(III. 20) 


where 


KA 
C_ 


cai 


cm  — =-  -*■ 


sec  cm   °C 


cal    sec 


T   ' 

2 


2m  c   gm   °C    cal 


C    sec  cal  gm  °C    sec 


2 .   Heat  Exchanger 

Assuming  no  heat  loss  to  the  outisde  the  following 
relation  holds:  The  total  change  in  heat  of  the  secondary 
sodium  within  the  heat  exchanger. 


AQ1  =  AQ2 


(III. 21) 


where    AQ   =  m  c(TBIl  -  TB01) 

=  the  heat  transferred  from  the  primary  sodium 
AQ2  =  KA  (TBI  -  TB2) 

=  the  heat  transferred  to  the  secondary  s 
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Using  average  temperature  of  primary  and 
secondary  side. 

Equation  III. 21  becomes  now 

KA(TB1  -  TB2)  =  m  C (TBI1  -  TB01)  (III. 22) 

or,  after  using  the  relation 


TBI  =  TBI1+TB01 


TBI1-TB2  =  ~    (TBI1-TB1) 

Ml 


.        ...    2  m  c 
where    Ml  = 


2 

gm    cal   sec  cm   °C   1 


KA     L sec   gm  °C      cal       2 

3  cm 


■*  dimensionless 


If  now,  as  done  before,  the  heat  flow  rate,  heat  capacitance 
of  the  primary  sodium  and  the  overall  transfer  coefficient  is 
assumed  constant,  Ml  may  be  seen  as  a  heat  exchanger  constant. 
The  dynamic  equation  for  the  secondary  side  is  derived 
as  in  section  II.B.l.b. 

n     r      dTB2  +  n 
Qp  -  CS  -dT"  +  QS 

where    Qp  =  KA  (TBI  -  TB2) 

Qs  =  m  c  (TH2  -  TC2) 

CQ  =  heat  capacitance  of  secondary  sodium 


or  after  rearranging 


dTB2  =  -   L_  +  -ip  TB2  +  —  TBI  +  -ip  TC2  (III.  23) 

dt        t3    13^        x3        x3 
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where 


KA 
C 


1 

sec 


1       2  m  c     1 

T3'  -      Cg     |_SeC_ 

3 .   Steam  Generator 

Since  the  flow  of  secondary  sodium  through  the  steam 
generator  can  be  changed  by  the  bypass  valve  in  the  secondary 
circuit  the  flow  rate  of  sodium  through  the  steam  generator 
will  be  defined  as  x*m,  where  x  is  the  valve  opening 
(0  <  x  <  1)  . 

x  =  1   valve  fully  open  (100%  flow  rate) 

x  =  0   valve  completely  closed 
As  in  section  II. B. 2.  the  following  equation  is  obtained 


KA(TB3  -  TS)  =  x  m  c(TBI2  -  TB02) 


eliminating  TB02  and  rearranging  gives 


(III. 24) 


TB3  -  TS  =  £-  (TBI2  -  TB3) 
M2 


where 


M2  = 


2  m  c 

KA 


=  a  dimensionless  steam  generator  constant 
under  the  assumptions  of  section  II. B. 2. 
For  the  dynamic  equation  of  the  steamside,  the  principle  used 
in  section  II.B.l.b.  is  again  applied  and  the  following 
equation  obtained 

dT^ 


QP  =  c 


s  d—  +  Qs 


(III. 25) 
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where    Q   =  KA  (TB3  -  TS) 

Cq  =  heat  capacitance  of  the  steam 

Qq  =  heat  consumed  by  the  steam  plant 
From  the  Carnot  cycle  it  is  known ,  that  the  work  done  in  a 
thermodynamic  machine  is  proportional  to  the  difference  of 
the  added  and  rejected  heat. 

H  a  (QA  -  QR) 
but      Qs  =  QA  -  QR 

and  Q   is  proportional  to  the  entrance  temperature  to  the 

steam  plant  and  Q   is  proportional  to  the  steam  outlet  tem- 
pi 

perature.   Taking  the  average  temperature  it  is  seen  that 

Q   ex  TS  Ci  W 

from  this  relation  Q   may  be  written  as  . 

Q   =  K  •  ty    •  TS 
K  =  constant,  assuming,  constant  flowrate,  heat 

capacitance  and  linear  temperature  change  with 
respect  to  pressure  over  the  range  of  operation 
cal 


°C 
ty   =   dimensionless  proportionality  factor  for  the 

load  demand  (0  <  ij;  <  1) 
y   =  i  ■+  100%  load 
ip   =  0  ->  0  load 
equation  (III. 25)  is  then 


KACTB3  -  TS)  =  Cs  |~  +  K  •  $    •  TS 
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and  after  rearranging 

dTi  1  1 

|i£.  =._  (_£  +  k'tJj)   TS  +  i-  TB3  (III. 26) 

5  5 

where    K1  =  ■=— 


Peal   °C    ,.     . 

0„  —v-  -*■  dimensionless 


i_  =  ka  r  i  ] 

t5    Cs  [secj 


C.   HEAT  TRANSPORT 

To  obtain  dynamic  equations  which  include  the  time  delays 
present  within  the  heat  transport  system  of  the  plant,  the 
Taylor  series  expansion  will  be  used.   Throughout  .this  section 
it  will  be  assumed,  that  there  is  no  heat  loss  during  the  heat 
transport  from  one  plant  component  to  another. 

1 .   Primary  Circuit 

The  transport  delay,  i.e.,  the  time  it  takes  for  one 
unit  volume  element  of  sodium  to  get  from  the  exit  of  the 
reactor  (with  temperature  TH1)  to  the  entrance  of  the  heat 
exchanger  (still  with  temperature  TH1,  but  now  called  TBI1)  is 
designated  T4.   Thus 

TH1  =  TBI1 
but   THl(t)  =  TBI1  (t  +  T4)  (III. 27) 

Expanding  the  right  side  of  equation  (III. 27)  in  a  Taylor 
series 

TBIl(t+T4)  =  TBIl(t)  +  T4  dTB^(t)  +  H.O.T.         (III. 28) 
where  H.O.T.  =  higher  order  terms. 
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Assuming  slow  temperature  changes,  which  is  reasonable 
for  a  large  stationary  power  plant,  all  higher  order  terms  may 
be  neglected. 

Replacing  TH1  and  rearranging,  equation  (III. 27)  may 
be  written  as 

2±|±±  =  -  ~  TBI1  +  ±-  TAV  -  ±-  TCI  (III. 29) 

For  the  transport  delay  from  the  heat  exchanger  back  to  the 
reactor  the  same  time  duration  will  be  assumed. 
Therefore   TB01  =  TCI 

but   TB01(t)  =  TCI  (t  +  T4)  (III. 30) 

Using  again  Taylor  series  expansion,  replacing  TBOl  and 
rearranging,  the  dynamic  equation  for  the  heat  transport  from 
the  heat  exchanger  back  to  the  reactor  is  given  as 

r)TT1       12  1 

2±^-  =  -  i-  TCI  +  £-  TBI  -  ±-  TBI1  (III. 31) 

at-       i4        14         4 

2  .   Secondary  Loop 

The  heat  flow  diagram  for  the  secondary  sodium  loop  is 
shown  in  Figure  11  where  Q   (n  =  1,2,... ,5)  is  the  rate  at 
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Fig.  11.  Heat  flow  diagram  for  secondary 
circuit. 


which  heat  is  transported  by  a  unit  volume  element.   The 
function  of  the  bypass  valve  Z  is  to  increase  or  decrease  the 
amount  of  sodium  flow  through  the  steam  generator:  it  does 
not  influence  the  temperatures.   This  may  be  shown  as  follows 

using  the  equivalent  of  an  electric  circuit,  the  "model 
equation"  for  node  A  is  written: 


Qo  =  Qn  -  Q- 


(III. 32) 


X 


m   c   TBI2  =  m   c   TH2  -  (1  -  x)m   c   TH2 


x  TBI2  =  TH2  -  (1-x)  TH2 
TBI2  =  TH2 
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Therefore  the  same  procedure  as  in  section  III.C.l,  may  be 
used  to  arrive  at  a  dynamic  equation  for  the  heat  transport 
delay  from  the  heat  exchanger  to  the  steam  generator.   The 
time  duration  of  the  transport  delay  will  be  called  T^ . 

Since    TH2     =  TBI2 

and     TH2(t)  =  TBI2  (t  +  Tc) 

and,  again  using  Taylor  series  expansion  for  TBI2 (t  +  T^)  and 
replacing  TH2  (t)  gives 

dTBI2  (t)  =  _  1_  TBI2  +  2_  TB2  _  1_  TC2         (Hi. 33) 
5  5         5 

For  the  heat  transport  from  steam  generator  to  heat  exchange 
the  "nodel  equation"  for  "node"  B  is 

Q5  =  Q3  +  ^4  UAA.J4J 

m  c   TC2  =  x  m   c   TB02  +  (1-x)  m   c   TH2 

Assuming  that  the  duration  of  the  transport  delay  is  again  T5 
for  the  flow  through  the  bypass  and  for  the  flow  from  the 
steam  generator,  equation  (III. 34)  becomes 

TC2(t  +  T5)  =  x  TB02(t)  +  (1-x)  TH2 (t)  (III.  35) 

Replacing  TB02  and  TH2 ,  using  the  relations 


M 
m2 

TBI2  -  TB3 

X 

TB3  -  TS 

'B3 

TBI2  +  TB02 

2 

'T3  0 

TC2  +  TH2 
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and  expanding  TC2 (t  +  T  )  in  a  Taylor  series,  equation  (III. 35) 
becomes 


dTC2      1        '  1  1        1 

S^-  =  ~  7TT   TC2  +  t^V  TB2  +  —±~   TS  +  — (7  TBI2    (HI.  36) 
at       A5  5  5  5 


where 

5       2-x 

T5 


5       2-2x 

T.  (x+M2) 


T         5 


5        2x  M2 


T  (x+M2) 


5       x(x-M2) 
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IV.   SYSTEM  SIMULATION 

A.   ANALOG  COMPUTER  PROGRAMMING 

The  coefficients  of  the  differential  equations  developed 
in  section  II  were  calculated  for  a  steady  state  solution  (see 
Appendix  C) .   Using  the  numerical  values  of  these  coefficients 
the  following  set  of  equations  which  describe  the  dynamic 
behaviour  of  the  plant  was  obtained: 
D  -         .58(N-D) 
TF   =  -  .5  TF  +  .5  TAV  +  725  N 
TAV   =  -  .5333  TAV  +  . 5  TF  +  4.8333  TCI 
TCI   =  -  .25  TCI  -  .0988  TBIl  +  .3488  TB2 
TBIl   =  -  .25  TBIl  +  .5  TAV  -  .25  TCI 

TB2   =  -6.9626  TB2  +  3.0233  TBIl  +  3.9394  TC2 
TBI2   =  -  .25  TBI2  +  .5  TB2  -  .25  TC2 
TC2   =  -  .2625  TC2  +  .275  TS  -  .0375  TB2  +  .025  TB2 
TS   =  -  .0421  TS  -  .0168  ■ ty    •  TS  +  .0421  TBI2 
N   =   D  /  (1-K.j,) 
To  simulate  these  equations  on  an  analog  computer  and  to  obtain 
maximum  accuracy,  the  gain  coefficients  must  to  be  scaled  with 
respect  to  the  maximum  voltage  level  of  the  analog  computer 
(see  Appendix  D) .   The  scaled  equation  used  for  the  simulation 
are  given  below. 

D   =    .58  CN-D) 
TF   =  -  .5  TF  +  .1833  TAV  +  .4  833  N 
TAV   =  -5.3333  TAV  +  1.3636  TV  +  4.3939  TCI 
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TCI   =  -  .25  TCI  -  .1235  TBI1  +  .3314  TB2 
TB*I1   =  -  .25  TBI1  +  .44  TAV  -  .2  TCI 
TC2   =  -  .2625  TC2  +  .2895  TS  -  .0474  TBI2  +  .025  TB2 
TB2   =  -6.9626  TB2  +  3.  +  78  TBI1  +  3.9394  TC2 
TBI 2   =  -  .25  TBI2  +  .3958  TB2  -  .1979  TC2 

TS   =  -  .0421  TS  -  .0168  •  4  •  TS  +  .0505  TBI2 
N   =  D/Cl-I^) 
The  temperature  variables  for  example  are  no  longer  in  °F  but 
in  volts . 

A  block  diagram,  using  both  the  symbols  shown  below  and  a 
linear  flowgraph  for  representing  the  analog  computer  setup 
may  now  be  drawn  (Figure  13,  14).   (See  Appendix  E  for  the 
analog  computer  diagram  corresponding  to  the  flowgraphs  shown 
in  Figure  14) 
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Figure  12.   Symbols  for  analog  computer  components 

It  proved  very  valuable  to  write  a  simple  Fortran  executive 
program  which  performed  the  above  calculations  each  time  a 


37 


kfA    kr 


FigD  13.    Analog  Computer  setup  -  reactor  core 
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test  or  simulation  run  was  made.   This  program  contained  the 
following  steps: 

a)  Read  in  steady  state  plant  variable  values, 
constants  and  chosen  maximum  plant  variable 
values 

b)  Calculate  the  numerical  value  for  the  coefficients 
of  the  differential  equations 

c)  Calculate  magnitude  scaling  factors 

d)  Calculate  scaled  coefficients  for  the  analog 
computer  components 

e)  Calculate  initial  conditions  for  analog 
integrators 

f)  Write  out  all  calculated  values  and  the  scaled 
and  unsealed  differential  equations 

g)  Calculate  the  initial  steady-state  output  values 
of  the  amplifiers 

h)   Set  all  potentiometers 

i)   Scan  all  amplifiers  and  compare  with  calculated 
values ,  write  out  errors 
The  amount  of  time  saved  by  using  this  executive  program 
during  the  many  runs  performed  was  appreciable.   The  program 
had  the  following  advantages : 

a)  Equipment  failures  were  detected  immediately  and 
isolated  to  only  a  few  components  by  the  test  part  of 
the  program. 

b)  The  scaling  procedure  was  done  by  the  computer.   Since 
the  maximum  plant  variable  values  were  not  known,  the 
optimum  scaling  coefficients  had  to  be  found  by  trial 
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and  error  in  many  runs.    By  just  changing  one  statement 
in  the  program  Cmaximum  variable  value)  the  analog 
computer  was  ready  for  the  next  run  within  seconds. 
c)   The  documentation  was  done  automatically  by  the  digital 
computer. 

B.   DYNAMIC  CHECK 

A  theoretical  curve  describing  the  behaviour  of  the  reactor 
of  the  proposed  plant  following  an  accidental  reactivity  input 
was  available  from  Ref.  [2].   The  validity  of  the  dynamic 
equations  describing  the  reactor  was  demonstrated  with  the  aid 
of  this  curve. 

To  check  the  solution  of  the  analog  computer  the  system 
equations  were  also  solved  on  a  digital  computer  (IBM  360/67) 
using  a  fourth  order  Runge  Kutta  routine  with  accuracy  testing 
of  single  and  double  increment  and  double  precision  calculation. 
The  results  of  these  dynamic  checks  are  stated  below. 

1.   Controlled  Accidental  Reactivity  Input 

The  accident  considered  assumes  that  two  control  rods 
drop  into  the  core.   This  would  correspond  to  a  reactivity 
insertion  of  100  $/sec.   It  is  further  assumed  that  scram 
Action  takes  place  .2  sec.  after  the  accident  starts. 

A  solution  of  the  system  equation  to  this  accident  was 
obtained  on  the  digital  computer  with  the  above  described 
routine . 


Because  of  the  nonlinearities  in  the  plant  dynamic  equa- 
tions, an  analytical  solution  for  the  calculation  of  the 
maximum  values  cannot  readily  be  obtained. 
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The  digital  computer  solution  gave  exactly  the  same 
results  for  the  temperature  variables  TAV  and  TF.   The  numerical 
values  for  the  neutron  density  N(t)  varied  slightly  in  magni- 
tude and  in  the  time  of  occurance.   This  slight  difference  can 
be  explained  by  the  simplifying  assumptions  which  were  made  in 
deriving  equation  II. 5  and  II. 6.   The  one  group  delayed 
neutron  and  the  prompt  jump  approximation  effects  the  time  of 
occurance  of  the  maximum  value  of  N(t).   By  neglecting  small 
reactivity  coefficient  (equation  11.10)  the  magnitude  of  N(t) 
is  influenced.   In  spite  of  this  small  difference  in  the 
result,  the  equations  used  for  the  simulation  of  the  reactor 
core  are  thought  to  be  accurate,  especially  since,  as  stated 
above,  the  temperature  TAV  and  TF  followed  very  precisely  the 
given  data. 

2 .  Change  in  Load 

The  first  part  of  the  check  on  the  accuracy  of  the 
analog  computer  solution  was  made  by  obtaining  a  digital  compu- 
ter solution  to  a  load  change  in  the  form  of  a  step  change  of 
magnitude  -  20%  for  20  sec.   The  results  are  given  in  Table  I 
and  in  Figure  15 . 

3 .  Change  in  Reactivity 

A  second  digital  computer  check  of  the  analog  computer 
solution  consisted  in  introducing  a  disturbance  in  the  reac- 
tivity input  in  the  form  of  a  +20  percent  reactivity  step 
change  for  20  sec.   The  results  are  given  in  Table  II  and  in 
Figure  16 . 
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Time  (sec)  TCIC°F)   TB2  [°F)    TS  (°F) 


0 

800.2 

885.1 

750.2 

2 

801.1 

886.6 

762.3 

4 

803.9 

890.4 

773.1 

6 

808.2 

895.5 

782.9 

8 

813.3 

901.5 

791.9 

10 

820.2 

908.6 

800.2 

15 

835.3 

924.9 

818.5 

20 

840.1 

941.5 

834.7 

22 

844.6 

946.5 

827.6 

24 

852.2 

948.5 

821.7 

26 

857.4 

948.8 

816.9 

28 

858.3 

948.1 

812.9 

30 

860.6 

946.8 

809.2 

35 

855.2 

942.2 

804.0 

40 

850.4 

937.3 

800.0 

TABLE  I 

Digital  computer  results;  temperature  changes  due  to 
step  change  in  load  of  -20%  for  0  <  t  <  20  sec. 
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Fig.   15.    Analog  computer  results  Load  change  as  in  table  I. 
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TimeCsec)  TCI  t°F)    TB2  C°F)    TSC°F) 


0 

800.1 

885.1 

750.2 

2 

800.3 

887.8 

750.2 

4 

801.1 

891.6 

750.5 

6 

802.2 

893.9 

751.0 

8 

803.3 

895.5 

751.8 

10 

804.3 

896.8 

752.8 

15 

806.5 

899.4 

755.8 

20 

808.9 

901.8 

758.9 

22 

809.7 

900.6 

760.2 

24 

810.0 

897.9 

761.0 

26 

810.0 

896.4 

761.7 

28 

810.0 

895.6 

762.0 

30 

810.0 

895,  3 

762,0 

35 

810.4 

895.2 

761.4 

40 

810.6 

895.0 

760.5 

TABLE 

II 

Digital  computer  results;  temperature  changes  clue 
to  step  change  in  reactivity  input  of  +  20  cents 
for  0  <  t  <  20  sec. 
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Fig.   16.    Analog  computer  results  Reactivity  change  as 
in  table  II. 
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4 .   Discussion,  of  the  Dynamic  Check  Results 

It  can  be  observed  from  the  above  results  that  the 
developed  equations  are  valid  for  the  plant  under  consideration, 
Comparison  of  the  given  digital  computer  solution  and  the 
analog  computer  graphs  shows  that  the  analog  computer  gives  a 
solution  which  is  within  the  bound  of  less  then  one  percent 
exact  (+  . 5%) . 

The  results  of  the  tests  have  demonstrated  the  validity 
of  the  simulation  within  the  above  bound  of  accuracy  and  under 
the  earlier  made  assumptions. 
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VI.   CONTROLLER 


A.   CONTROL  STRATEGY 

The  inherent  incapability  of  the  plant  to  follow  load 
changes  which  is  due  to  the  small  reactivity  feedback  coef- 
ficient of  the  coolant,  is  a  disadvantage  in  the  control  of  a 
Fast  Breeder  Reactor.   It  makes  complete  control  of  the  coolant 
temperature  in  the  reactor  necessary.   An  automatic  controller 
for  power  level  control  must  be  developed  whose  two  main 
purposes  are 

a)  coolant  temperature  level  must  be  kept  constant 

to  obtain  optimum  breeding  ratio  and  optimum  perfor- 
mance of  the  reactor  and 

b)  steam  conditions  are  to  be  kept  constant  for  optimum 
performance  of  the  steam  plant. 

On  the  other  hand  several  advantages  which  are  inherent 
for  the  plant  help  the  controlling  action  to  achieve  smooth 
and  efficient  power  level  changes.   The  reactor  is  automati- 
cally isolated  from  the  load  by  the  secondary  sodium  loop  and 
to  a  certain  degree  by  the  primary  sodium  loop.   Therefore, 
the  inlet  temperature  to  the  core  changes  only  very  slowly. 
These  two  circuits  also  serve  as  a  heat  reservoir.   This  allows 
a  temporarily  larger  rate  of  load  change  than  the  rate  of  power 
change  which  the  reactor  can  achieve. 

The  rate  of  power  change  obtained  from  the  reactor  is 
limited  by  the  maximum  rate  of  change  of  reactivity.   This 
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cent  s 
value  is  taken  as  2     /sec.  [2]   It  is  also  limited  by 

mechanical  factors  such,  as  thermal  stress  within  the  reactor, 
heat  exchanger  or  steam  generator  and  speed  of  mechanically 
driven  valves.   The  thermal  stresses  were  not  investigated 
here,  it  can,  however,  be  stated  that  in  none  of  the  test  runs 
the  limiting  values  for  the  rate  of  change  of  temperature, 
e.g.  the  coolant  in  the  core  with  100°F/sec. ,  as  given  in 
Ref .  [2] ,  was  observed.   The  limiting  factor  in  this  investi- 
gation was  found  to  be  the  rate  of  change  of  reactivity. 

The  control  strategy  is  based  on  the  following  thoughts : 
Because  of  the  extreme  long  time  constant,  a  pure  pro- 
portional controller,  which  acts  only  following  the  occurance 
of  errors  in  the  to-be-controlled  quantities  was  not  sufficient. 
Therefore,-  an  anticipatory  or  derivative  control  had  to  be 
implemented.   To  get  the  best  derivative  control  action  signal, 
the  derivative  or  the  rate  of  change  of  temperatures  which  had 
to  lead  to  changes  in  the  to-be-controlled  temperature  had  to 
be  sensed  at  the  earliest  possible  time.   This  steady  state  - 
or  at  most  very  slow  power  level  change  -  strategy  is  not 
valid  for  larger  power  level  changes,  in  this  case  preprogrammed 
controlling  action  had  to  be  implemented. 

B.   CONTROLLER  DESIGN 

The  choice  was  given  between  a  pure  analog-logic  controller 
or  a  digital  controlling  program. 

The  analog  controller  had  the  advantage  that  it  could  have 
been  implemented  comparitively  easy  for  a  plant  simulation  of 
this  size.   On  the  other  hand  it  was  felt  to  be  more  re     tic 
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to  design  an  automatic  controller  which,  uses  the  digital 
computer  since  modern  power  plants  are  equipped  with  digital 
computers  used  for  many  purposes  -  from  bookkeeping  of  the 
plant  datas  to  paycheck  writing  -  that  a  digital  controller 
would  represent  only  a  small  and  comparitively  low  cost 
addition  to  the  main  computer.   Investigations  have  also  shown 
[8],  that  for  large  scale  plant  controllers  digital  controllers 
become  more  economical  than  analog  controllers.   Therefore,  it 
was  decided  to  design  a  digital  controller.   The  following 
controlling  devices  were  "built"  in  to  the  plant: 

a)  control  rods.   It  was  assumed  that  the  rate  of  change 
of  reactivity  could  be  varied.   For  the  preprogrammed 
controlling  action  the  maximum  rate  of  2  cents/sec.  was 
assumed;  for  the  steady  state  controlling  action  a 
value  of  .1  times  the  maximum  rate  was  found  by 
experiment  to  give  the  smoothest  control  dynamics . 

b)  Bypass  valve  in  secondary  circuit.   The  opening  time 
(with  maximum  speed)  for  the  preprogrammed  phase  was 
assumed  to  be  20  sec.  from  completely  closed  to  fully 
open;  for  the  small  controlling  action  .2  of  this  value 
was  used. 

Because  of  the  constantly  changing  value  of  the  bypass  valve 
opening,  several  gain  factors  in  the  dynamic  equations 
developed  in  section  II  changed  during  operation.   As  a  result 
of  this  the  analog  computer  setup  had  to  be  modified  (see 
Appendix  E) . 
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The  above  described  controlling  devices  where  controlled 
by  the  digital  computer  which  in  turn  received  the  data  from 
the  simulated  plant  on  the 


Control 
Rods 


~/*~\   TAV 

(\  r  V  >-i 
-*1 


TCI 


l  I 


i  i 


IU-  - 


I 


Fig.   17.     Plant  and  control  points. 


analog  computer  (see  Fig.  17)  . 

The  control  strategy  called  for 

a)   preprogrammed  power  level  change  control.   The  control 
trajectory  for  this  control  was  obtained  by  introducing 
load  changes  and  correcting  the  errors  in  TCI  and  TS 
with  manually  operated  potentiometer,  which  simulat 
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control  rods  and  bypass  valve  opening.   The  resulting 
datas  were  evaluated  with  the  aid  of  a  least  square 
straight  line  approximation  program  to  give  the 
control  trajectory  for  the  control  reactivity  and  valve 
opening  as  a  function  of  the  load  (see  Appendix  F) . 

b)  Derivative  control. 

Following  the  thought  described  in  IV. 1,  the  rate  of 
change  in  TAV  and  TB2  were  calculated  by  storing  values 
every  10    of  a  second  and  then  calculating  the  rate  of 
change  over  one  second. 

c)  Proportional  control. 

The  to-be-controlled  temperatures  TCI  and  TS  were 
compared  every  10    of  a  second  with  the  set  value 
and  the  error  was  calculated . 
To  take  the  addition  of  noise  in  actual  plant  measured 

variables  into  consideration,  all  measured  values  were  put  into 

a  dead  zone.   This  had  also  the  advantage  of  smoothing  the 

controlling  action. 

A  flow  chart  of  the  basic  program  for  the  digital  controller 

is  given  in  Figure  18.   The  individual  parts  of  the  program 

are  described  in  detail  in  Appendix  G. 
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VI.   RESULTS  AND  COMMENTS 

The  results  which  were  achieved  by  controlling  the  plant 
with  the  described  controller  are  shown  in  this  section  as 
stripchart  recordings  in  Figure  19  through  Figure  21.   The 
test  inputs  were  step  changes  in  load  of  different  sizes  and 
at  different  time  intervals. 

Figure  19  shows  the  plant  dynamics  following  a  step  change 
in  load  of  -50%.   This  step  change  does  not  seem  to  be  very 
realistic,  it  proves,  however,  that  the  initial  overshoot  in 
the  to-be-controlled  temperatures  of  little  more  than  100 °F 
is  damped  down  very  effectively  to  a  secondary  overshoot  -  or 
first  undershoot  -  of  about  20°F.   The  initial  overshoot 
cannot  be  avoided  since  the  controlling  devices  work  at 
maximum  allowable  speed  during  the  preprogrammed  part.   The 
same  statements  hold  for  the  step  change  in  load  of  +50%  which 
was  introduced  after  the  transients  had  almost  settled.   Even 
a  load  change  of  this  amount  does  not  endanger  the  plant  -  due 
to  the  fast  working  and  direct  preprogrammed  controller  -  which 
may  be  seen  from  the  value  of  N  which  increases  only  about  12%. 

It  is  believed,  that  this  test  input  and  the  resulting 
transient  demonstrate  the  validity  and  necessity  of  the 
division  of  the  control  program  into  its  two  main  parts,  it 
also  demonstrates  effectiveness  of  the  preprogrammed  controller 

Figure  20  shows  the  plant  dynamics  following  a  step 
change  of  different  magnitudes  and  at  time  intervals  which 
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are  shorter  than  the  settling  time  of  the  dynamics.   A  20% 
step  change  in  load  produces  only  a  25° F,  or  35 °F  respectively, 
overshoot  in  the  controlled  temperature.   The  increase  in  N 
is  only  4%.   The  interplay  of  preprogrammed  and  small  deviation 
control  may  also  be  seen  from  Figure  20  by  considering  the 
reactivity  and  bypass  valve  curve. 

Figure  21  was  taken  over  a  time  of  more  than  1  hour. 
Though  step  changes  in  load  of  10%  do  not  seem  to  be  very 
practical  this  recording  shows  a  more  realistic  picture  of  the 
plant  controller.   It  especially  proves  the  accuracy  of  the 
steady  state  controller.   After  the  settling  of  the  dynamics 
this  controller  is  able  to  keep  the  error  in  the  controlled 
temperatures  within  a  bound  of  less  than  +  3°F  at  a  controlled 
temperature  level  of  800°F  for  TCI  and  750°F  for  TS . 

The  above  results  demonstrate  the  effectiveness  of  the 
designed  controller  for  the  plant  layout  simulated  in  this 
investigation.   Since  conventional  stability  calculations  like 
root-locus  plot,  Bode  plot  or  Nyquist  diagram  cannot  be 
applied  directly  to  a  controller  of  the  kind  developed  here, 
this  topic  is  not  included  in  this  report.   It  can  however  be 
stated  that  no  instabilities  following  any  of  the  applied  load 
changes  were  observed  during  the  simulation. 

The  addition  of  more  plant  components  in  the  simulation 
of  the  system  will,  as  stated  in  section  II,  not  change  the 
basic  form  of  the  dynamic  system  equation.   It  will  however, 
change  the  dynamic  response  to  a  certain  degree.   The  same 
is  true  for  the  nonlinearities  which  were  not  taken  into  account 
in  the  simulation  in  this  report. 
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Further  investigations  can  include  the  above  considera- 
tions, however,  a  simulation  of  a  complete  plant  which  includes 
all  components  and  takes  all  nonlinearities  into  account  is 
only  possible  on  a  very  large  scale  digital,  or  better,  hybrid 
computer.   If  a  computer  of  this  size  is  not  available  an 
accurate  simulation  can  be  done  component  by  component  as  in 
Ref.  [3]  for  a  LMFBR.   Accurate  stability  investigation  may 
then  be  performed  by  simulation. 
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APPENDIX  A 

A.   CALCULATION  OF  THE  REACTIVITY  FEEDBACK  COEFFICIENTS 
The  reactivity  coefficients  considered  are 

a)  Doppler  coefficient 

b)  Fuel  expansion  coefficient 

a)   The  Doppler  coefficient  is  a  nonlinear  function 
of  the  fuel  temperature  and  is  shown  in  Fig.  22. 


1000 


4000 


Fig.  22.  Doppler  coefficient. 


Since  the  operating  temperature  of  the  reactor  is  at  2400 °F  an 

-5 
average  doppler  coefficient  of  -.35*10    could  be  taken.   It 

can  be  seen,  that  even  at  temperature  variations  of  +400° , 
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which  have  to  be  avoided  in  a  fast  reactor,  the  change  in 

.  .     .  -5 

doppler  coefficient  is  at  most  +.5*10   .   To  account  for  a 

certain  degrees  of  uncertainty  a  conservative  numerical  value 

-5 
of  -.3.10    will  be  used.   Changing  this  value  to  dollar  units 

-4 
of  reactivity  with  g  taken  as  29.88x10    the  doppler  coefficient 

is  obtained  as 


-5 
Kd  =  — ~-3xl° =  -1.04xl0~3   dollars 

29.88xl0~4 


b)   The  fuel  expansion  coefficient  is  considerably 
smaller  than  the  doppler  coefficient  and  since  its  reliability 
for  axial  fuels  is  highly  uncertain  the  conservative  numerical 
value  of  -1.0 -10   will  be  taken.   The  conversion  into  dollar 
units  reactivity  gives: 


K„YT3  =  — 1-'10 =  _  .34xl0"3  dollars 

A 

29.88x10 


The  total  feedback  reactivity  is  therefore  given  as 


KF  ~  KD  +  KEX 


=  -.00104  -  .00034 

=  -.00138 

=  -.0014   dollars 
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APPENDIX  B 


A.   SOLUTION  OF  THE  STEADY  STATE,  SECOND  ORDER  DIFFERENTIAL 
EQUATION  FOR  THE  RADIAL  HEAT  FLOW  IN  A  CYLINDRICAL  FUEL 
ELEMENT 

The  heat  balance  equation  is  given  as 


Q  2irr  dr  L  =  Qr+dr  -  Qr  (B.l) 

The  differential  heat  quantities 

dQ 

Q.^   =  Q   +dQ   =  Q   +  -5-^-  dr  (b.2) 

*r+dr    r     r    r   dr 


dT  dT 

Q   =  -  K_  A  —  =  -  K.  2TTr  L  ~  (B.3) 

r      f   dr      f       dr 


equation  (B.2)  becomes  now 

9 

dT  d""T    dT 

Q  ^    =   "    2tt  L  K   p-  -  2tt  L  K(r  ?—=r  +   ~)  dr         (B.4) 
r+dr  r  dr  -  2   dr 

dr 


or 


2 

d  T    df 
Q^  2-rrr  dr  L  =  -  2n  L  K(r  2-=-  +  —■)    dr  (B.5) 

F  ,  2   dr 

dr 


which,  after  concellation  and  rearranging,  may  be  written  as 


dr 


This  differential  equation  can  be  solved  by  standard  methods, 

the  solution  is 

2 
T  =  -  Qp  ~  +  C1  In   r  +  C2  (B.7) 
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The  boundary  conditions  are 


dT 
dr 


B   0 

r=0 


=  T 
r=0 


=  T 
r=R 


Applying  these  boundary  conditions  gives 

qfr2 

T   -  T   = 


c     s     4K 


Since 


QR  =  iTr2  L  QF  (B.8) 


it  follow  that 


Qs  =  ttR2  l  Qf  (B.9) 


or,    combining    (B.8)    and    (B.9) 


Qc    =    477    K    L(T       -    T    )  (B.10) 


Setting  2tt  R  L  =  A,   the  surface  of  the  fuel  element 

Tc  "  TS 
Qs  =  2  Kf  A   °  R  (B.ll) 

Equation   CB.ll)  holds  for  the  region  inside  the  fuel.   For 

the  cladding  region  K  ,  ,  would  have  to  be  used  instead  of  Kf 

and  for  the  region  in  the  coolant  K    .       Since  average 

coolant 
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temperature  are  used  and  in  addition  the  individual  thermal 
conductivities  are  considered  constant,  an  average  thermal 
conductivity  will  be  taken  and  CB.3)  may  be  written  as 


r3T 
Q   =  -  K  2ur   L  ~  (B.12) 

*r  dr 


where   r  =   the  distance  measured  from  the  center  of  the  fuel 
element. 

T   is  taken  as  TF,  the  average  fuel  temperature  and   Tq 
as  TAV,  the  average  coolant  temperature. 

Equation  CB.5)  becomes  then 


Q„  =  KA  (TF  -  TAV)  (B.13) 

t 
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APPENDIX  C 

A.   CALCULATION  OF  THE  GAIN  COEFFICIENTS  OF  THE  SYSTEMS 
DIFFERENTIAL  EQUATIONS 

The  gain  coefficients  in  the  individual  equations  may  be 

solved  for  by  using  the  numerical  values  given  for  a  specific 

plant.   This  would  require  a  certain  amount  of  engineering 

approximation,  since  exact  values  are  not  available  for  all 

quantities.   Instead  of  following  this  approach,  only  given 

steady  state  temperature  levels  for  100%  power  (ip  =  1)  at  a 

bypass  valve  opening  of  95%  (X  =  .95)  and  given  time  constants 

were  used.   These  values  are  19] : 

a)  Steady  state  temperatures 
TF    =   2400°F 

TH1  =  1100°F  =  TBI1 
TCI  =  800°F  =  TB01 
Th2  =  1050°F  =  TBI2 
TC2  =  720°F  =  TB02 
TS    =    750°F 

b)  Time  constants 

Tl  =  2  sec  =  T2 

T3  =  4  sec  =  T4  =  T6 

T5  =  10  sec 

T7  =  .1  sec 

Additional  constants  used  [3]  : 
X     =   .58 
B     =   .002988 
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Setting  these  constants  into  the  obtained  dynamic  equation 
these  equations  were  "forced"  into  an  accurate  steady  state 

solution  by  solving  for  the  unknown  gain  coefficient  in  the 

dT 
steady  state  form  of  the  differential  equation  (-=—  =  0)  . 

This  procedure  is  demonstrated  with  one  example: 

The  equation  for  heat  transfer  within  the  core  was  obtained 
as  : 


TF  =  -  =-i  TF  +  =i  TAV  +  A'N 
1        1 


which  becomes  in  its  steady  state  form 


0  =  -  |  (2400)  +  |  (950)  +  A1  (1) 


and  solving  for  A'  gives 

A1  =  750 

which  is  the  desired  coefficient  for  N.  Using  this  approach 
for  all  equations  individually  leads  to  the  system  equations 
given  in  section  II. 
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APPENDIX  D 

A.   SCALING  OF  THE  SYSTEM  EQUATIONS 

As  mentioned  in  section  III  the  actual  equations  must  be 
scaled  with  respect  to  the  maximum  voltage  levels  of  the 
analog  computer.   It  is  desirable  to  get  maximum  output  values 
of  the  equation  variables  close  to  the  computer  maximum 
voltage  level  to  maintain  maximum  accuracy  in  the  solution  of 
the  differential  equations.   To  be  able  to  perform  this  scaling 
in  an  optimum  way  an  exact  knowledge  of  the  maximum  values 
of  the  plant  variables  is  necessary.   If  these  variables  are 
not  known  exactly,  optimum  scaling  is  done  by  trial  and  error: 
a  maximum  value  is  assumed,  the  equations  are  scaled  for   a 
computer  run  and  the  output  signal  is  observed.   If  this 
output  signal  is  either  to  high  or  too  low,  the  equations  have 
to  be  rescaled. 

The  magnitude  scaling  procedure  uses  scaling  coefficients 
which  convert  given  plant  variables  with  their  individual 
dimensions  to  computer  variables  with  the  deminsion  volt. 

Taking  the  example  of  TAV  it  is  written 


TAV      =   SLTAV  •  TAV 
where  SLTAV  =   scaling  coefficient  for  the  plant  variable 

TAV  with  dimension 


°F 


V 


TAV  =   computer  variable  for  average  coolant  temper- 
ature in  core  with  dimension  (V) 
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The  analog  computer  used  here  (COMCOR  Ci  5000)  had  a 
maximum  voltage  level  of  +  100V.   Assuming  a  maximum  value  of 
1100°F  for  TAV  the  scaling  coefficient  is  obtained  as 


SLTAV  -  TAVX  -  110°  -  i l 
SLTAV  "  TAVX  "  "TOO  ~  -11 


F 


V 


Using  this  procedure  the  following  numerical  values  for  the 
individual  scaling  coefficient  proved  to  be  best  suited: 

SLD  =     .02 

SLTF  =  30. 

SLTAV  =  11. 

SLTC1  =  10. 

STBI1  =  12.5 

SLTB2  =  9.5 

STBI2  =  12. 

SLTC2  =  9.5 

SLTS  =  10 . 

The  scaling  procedure  for  a  single  dynamic  equation 
proceeds  as  follows : 

Taking  equation  (11.20) 

TAV  =  5.3333  TAV  +  .5        TF  +  4.8333  TCI 

Replacing  the  plant  variables  by  computer  variables  gives 


SLTAV  •  TAV  =  -  5. 3333* SLTAV 'TAV  +  .5 'SLTF  •  TF 
+  4.833  -    SLTC1  •  TCI 

Dividing  through  by  SLTAV  and  using  numerical  values  gives 


TAV  =  -  5.3333  TAV  +  1.3636     TF  +  4.3939     TCI 
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Using  this  procedure  for  all  equations  and  steady  state 
values  gives  the  set  of  scaled  equation  given  in  III  and  the 
corresponding  initial  conditions  for  the  integrators  of  the 
analog  computer.   These  initial  conditions  were  found  to  be 

TF  (0)   =   80.      (V) 
TAV  (0)   =   86.3636  (V) 


TCL  (0)   =   80       CV) 


TBI1  (0)   =   88       CV) 


TB2  (0)   =   93.1579  (V) 


TBI2  (0)   =   87.5     CV) 


TC2  (0)   =   75.7895  (V) 
TS  (0)   =   75.      CV) 

The  steady  state  control  rod  reactivity  required  is  found 
to  be  3.3C  dollar  which,  after  being  scaled ,  corresponds  to 
33.6  (V) . 
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APPENDIX  E 


A.   ANALOG  COMPUTER  SETUP  FOR  CONTROLLER  IMPLEMENTATION 

It  was  shown  in  II,  that  several  gain  factors  in  the  system 
dynamic  equation  were  a  function  of  the  bypass  valve  opening 
X.   Since  the  bypass  is  used  as  one  of  the  controlling  ele- 
ments in  the  plant  and  therefore  continually  changes  its 
numerical  value,  a  way  had  to  be  found  to  continually  change 
the  gain  factors  which  are  functions  of  X.   For  the  simulation 
without  controller  the  gain  factors  were  implemented  as  com- 
puter set  potentiometers,  which  take  a  considerable  amount  of 
.time  to  be  set  or  adjusted.   A  continuous  dynamic  solution  or 
simulation  of  the  plant  was  therefore  not  possible  with  this 
setup. 

The  following  method  was  used  to  get  continuously  changing 
gain  factors  and  therefore  a  continuous  simulation: 

The  potentiometers  were  replaced  by  electronic  multipliers 
which  had  as  input  the  to  be   multiplied  variable  (temperature) 
and  the  output  of  a  trunkline  which  was  controlled  by  the 
digital  computer.   Figure  2  3 

tf>-  TS     \\/    K  -  j>  -TS^ 


A 


TRUNKLINE 
Set  to  K   [v] 


+100v 


Fig.  23.    Implementation  of  varying  gain 
factors  . 
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shows  the  analog  implementation  for  the  second  gain  factor 
in  the  differential  equation  for  TS  (K'^'TS) .   The  voltage 
representing  TS  is  multiplied  by  a  voltage  representing  ij> 
(0[V]  <  i>    <  100 IV]  corresponds  to   0  <  ty    <    100%  load).  ty 
was  changed  by  a  manually  set  potentiometer  on  the  analog 
board.   The  output  of  this  multiplier  is  fed  into  another 
multiplier  with  the  trunkline  as  the  second  input.   Each 
time  (scan  intervall  =  .1  sec)  the  controller  changed  X,  it 
simultaneously  calculated  all  dependent  gain  factors  and 
adjusted  the  corresponding  trunkline.   Due  to  the  electronic 
multiplication  continuous  adjustment  of  the  gain  factors  with- 
out any  timelag  was  thereby  achieved. 

Figure  24  shows  the  resulting  analog  computer  setup  for 
the  plant  part  without  the  rs?.ctor.   For  more  clarity  of  the 
graph  no  numerical  values  were  set  in  and  all  initial  condi- 
tion inputs  to  the  integrators  with  the  corresponding 
potentiometers  were  left  out.   Figure  2  4  corresponds  to  the 
linear  flowgraph  of  the  plant  system  equation  given  in  III. 

Since  there  are  no  changing  gain  factors  in  the  reactor 
part  of  the  simulation  Figure  13  in  III  is  not  changed. 
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•    Fig.  24.    Analog  computer  setup  for  plant, 
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APPENDIX  F 

A.   CONTROL  TRAJECTORY  FOR  PREPROGRAMMED  POWER  LEVEL 
CHANGE  CONTROL 

300  data  points  were  obtained  by  changing  the  load  over 

the  whole  range  from  100%  to  0%  in  small  steps.   For  these 

runs  the  analog  computer  was  in  1000  :  1  timescale  mode, 

therefore  the  steady  state  error  for  a  certain  load  was 

immediately  available.   (see  Figure  27).   Using  hand  set 

potentiometers  voltages  "were  supplied  to 

a)  the  analog  computer  which  corresponded  to  control 
reactivity 

b)  the  digital  computer  which  corresponded  to  the 
bypass  valve  position. 

The  digital  computer  calculated  the  variable  gain  factors  and 
adjusted  the  trunlines. 

The  "control  reactivity  voltage"  and  the  "bypass  voltage" 
were  adjusted  until  the  steady  state  error  was  zero.   The 
obtained  data  points  were  then  used  to  obtain,  with  the  aid 
of  a  least  square  polonomial  fit  program,  straight  line 
approximations  as  precomputed  control  trajectory  for  load 
changes  in  the  plant.   These  data  points  are  shown  in  Figure 
25  and  Figure  26. 

As  may  be  seen  from  Figure  27,  it  is,  with  the  controlling 
aids  of  this  plant  simulation  possible  to  control  TCI 
completely  over  the  whole  load  range,  TS  can  only  be  controlled 
over  a  load  range  from  10  0%  load  to  about  20%  load. 
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state  error  for  TS  at  zero  load  is  about  +80'°F.   This  is  not 
a  restriction  on  the  controllability  of  the  plant. 
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J&^ 


x.  "*<  K  ££$$£  ^- 


010 


Fig.  25.    Experimental  data  points. 

Control  Trajectory  of  reactivity 
Scale    X:    1  inch  =  10%  load 

Y:    1  inch  =  dollar  reactivity 


c  Km-}* 


n 

K 


3.13 


Fig.  26.     Experimental  Data  points. 

Control  Trajectory  of  bypass  opening 
Scale:     X:    1  inch  =  10%  load 

Y:    1  inch  =  10%  valve  opening 
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APPENDIX  G 

A.   DIGITAL  CONTROL  PROGRAM 
1.   General 

The  scan  interval  was  chosen  as  .1  sec,  i.e.,  the 
digital  computer  went  through  the  designed  control  loop  every 
.1  sec.  and  performed  the  necessary  control  action  by  scanning 
the  assigned  variables,  performing  the  necessary  calculations 
to  get  the  derivatives  or  errors,  and  then  going  through  all 
the  logic  steps,  implemented  as  Fortran  IF-statements .   If 
control  action  was  necessary  it  was  performed  by  calling  the 
following  subroutines. 

a)  REACP :    Increase  reactivity  at  the  rate  of  2  cents/ 
sec.   Because  of  the  scan  interval  the  voltage  of  the 
trunkline  which  represented  reactivity  input  was 
increased  by  an  amount  which  corresponded  to  .2  cents/ 
.1  sec.   Inspite  of  being  a  step  change  the  short 
time  interval  in  which  these  changes  were  performed 
made  the  controlling  action  an  almost  continuous  one. 
The  time  delay  of  an  electric  motor  which  would  be 
used  to  perform  this  controlling  action  in  a  real 
plant,  was  neglected.   This  is  justifiable  in  view 

of  the  short  time  constant  of  a  motor  with  respect 
to  the  time  constants  of  the  plant. 

b)  REACM:    Decrease  reactivity  at  a  rate  of  2  cents/sec. 
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APPENDIX  G 

A.   DIGITAL  CONTROL  PROGRAM 
1.   General 

The  scan  interval  was  chosen  as  .1  sec,  i.e.,  the 
digital  computer  went  through  the  designed  control  loop  every 
.1  sec.  and  performed  the  necessary  control  action  by  scanning 
the  assigned  variables,  performing  the  necessary  calculations 
to  get  the  derivatives  or  errors,  and  then  going  through  all 
the  logic  steps,  implemented  as  Fortran  IF-statements .   If 
control  action  was  necessary  it  was  performed  by  calling  the 
following  subroutines . 

a)  REACP :    Increase  reactivity  at  the  rate  of  2  cents/ 
sec.   Because  of  the  scan  interval  the  voltage  of  the 
trunkline  which  represented  reactivity  input  was 
increased  by  an  amount  which  corresponded  to  .2  cents/ 
.1  sec.   Inspite  of  being  a  step  change  the  short 
time  interval  in  which  these  changes  were  performed 
made  the  controlling  action  an  almost  continuous  one. 
The  time  delay  of  an  electric  motor  which  would  be 
used  to  perform  this  controlling  action  in  a  real 
plant,  was  neglected.   This  is  justifiable  in  view 

of  the  short  time  constant  of  a  motor  with  respect 
to  the  time  constants  of  the  plant. 

b)  REACM:    Decrease  reactivity  at  a  rate  of  2  cents/sec. 


81 


c)  RACP  and  RACM:    same  as  under  a)  and  b)  but  with 
the  rate  divided  by  10. 

Subroutine  REACP  and  REACM  were  called  if  controlling 
action  was  performed  by  preprogrammed  controller. 
Subroutine  RACP  and  RACM  were  called  if  controlling 
action  was  performed  by  small  deviation  controller. 

d)  VENTP :    Open  bypass  valve  at  described  rate  (0  -  100% 
in  20  sec).   The  subroutine  changed  the  value  of  X, 
calculated  the  new  gain  factor  and  adjusted  the 
corresponding  trunklines  to  the  electronic  multipliers . 
The  same  statements  made  in   a)  with  respect  to 
continuous  action  and  time  constant  hold  also  here. 

e)  VENTM:    Decrease  bypass  valve  opening  at  the  described 
rate . 

f)  VETP  and  VETM :    same  as  under   d)  and   e)  but  with 
rate  divided  by  5. 

Subroutine  VENTP  and  VENTM  were  called  by  the  pre- 
programmed controller  subroutine  VETP  and  VETM  by  the 
small  deviation  controller. 
The  digital  controller  program  was  divided  in  3  parts: 

a)  Routing 

b)  Preprogrammed  controller 

c)  Small  deviation  or  steady  state  controller. 
These  3  parts  are  described  below. 

2 .   Routing 

The  decision  whether  to  go  into  the  preprogrammed 
part  of  the  controller  or  to  the  small  deviation  controller 
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was  based  on  the  rate  of  change  of  the  load.   For  a  load 
change  of  less  than  .2%/sec.  the  control  action  was  performed 
by  the  small  deviation  controller.   The  flowchart  of  the 
routing  part  is  given  in  Fig.  28. 
3  .   Preprogrammed  controller 

The  purpose  of  this  controller  was,  to  bring  the 
controlling  devices  directly  to  their  set  points,  i.e.  to  the 
points  which  were  calculated  by  the  straight  line  approximation 
-  as  a  function  of  the  load  -  to  lie  on  the  control  trajectory. 
This  controller  performed  its  function  without  regard  of  the 
dynamics  of  the  plant.   The  load  was  scanned,  the  set  point 
calculated  and  the  respective  subprogram  (Controlling  action) 
called.   Fig.  29  shows  a  flow  chart  of  this  part  of  the 
digital  controller. 

4 .   Small  Deviation  or  Steady  State  Controller 

As  stated  in  III  this  part  of  the  controller  used 
proportional  and  derivative  control.   The  decisions  of  the 
digital  computer  were  based  on  the  logic  tables  given  as 
Table  III  and  Table  IV.   The  procedure  used  here  is  the  same 
as  in  analog  logic  circuit,  instead  of  writing  expressions  in 
Bolean  algebra,  Fortran  IF-statements  were  written  for  the 
logic  decisions  of  the  digital  computer. 

Controlling  action  for  change  of  reactivity  was  based 
on  TAV  and  ETCl  (error  in  TCI) .   The  logic  table  used  is 
Table  III. 

Controlling  action  for  change  of  bypass  valve  value 
was  based  on  TB2  and  ETS .   The  loqic  table  used  is  Table  IV. 
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With  respect  to  earlier  thoughts  any  decision  made  on 
the  basis  of  Table  III  was  weighted  higher  and  if  controlling 
action  was  required,  Table  IV  was  neglected,  i.e.  the  controller 
tried  first  to  bring  ETCI  to  zero  and  keep  it  zero  before  ETS 
was  controlled.   This  was  done  because  of  the  dynamic  of  the 
plant.   No  logic  scheme  could  be  arrived  at,  which  would  have 
made  combined  control  possible.   It  will  be  shown  in  section 
V,  that  the  performance  of  this  controller  was  very  good.   A 
flow  diagram  of  this  part  of  the  digital  controller  is  given 
in  Figure  30 . 
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Start 


JL 


Array  Load  (I) 

I  =  1,10 
Scan  Freq .   10  cps 


X 


Load  (I)  =  Load  (1+1)  << 


1L 


Scan  Freq  .  10  cps 


Scan  Load  (10) 


A. 


A    Load  = 
Load  (10)  -  Load  (1) 


>  .2%/sec 


Preprogrammed 
Controller 


Small  Deviation 
Controller 


Fig.  28.    Routing  part  of  digital  controller. 
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Routing  Part 
of  Controller 


e 


Reactivity 
Bypass 


=  F(Load) 


Call  subroutine  for 
necessary  control 
action 


h 


_L 


Scan  Load 

x 


-^      Scan  Freq  .  =10  cps 


Fig.  29.     Preprogrammed  controller  part  of  digital  controller. 
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Logic  Tables  for  Small  Deviation  Controller 


TAV 
0 

+ 


0 

+ 

0 

+ 


TCI 

Reactiv 

0 

./. 

0 

D 

0 

I 

+ 

D 

+ 

D 

+ 

•A 

^ 

I 

- 

./• 

pm 

I 

Table  III.   Reactivity  control. 


Symbols : 


TB2 


0 

+ 

0 

+ 

0 

+ 


ETS 


B.  P.  valve 


0 
0 
0 

+ 
+ 
+ 


D 
I 
D 
D 

./• 

I 

./• 
I 


Table  IV.   Bypass  valve  control. 


+    Error  or  Deriv.  >    0 
Error  or  Deriv.  <  Q 

0  Error  or  Driv.   =  0 
D    Decrease 

1  Increase 
./.   No  action 
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Routing  or  Preprogram- 
ed Part  of  Controller 


Yes 


ETS(I)=ETS(I+1) 


\_ 


Call  Scan 
ETS 


I 


ETS(10)=ETS 


Array     I  =  1,10 
ETS  (I) 
ETC  1  (I) 
TAV(I) 
TB2(I) 


-^ 


Calculate  Derivative 
TAV,  TB2 


Logic  Tables 


To  Routing 
Part 


Action 


Action 


A 


Fig  „  30.     Small  deviation  controller  part  of 
digital  controller. 
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APPENDIX  H 

A.   DIGITAL  COMPUTER  PROGRAM 

The  following  pages  are  a  copy  of  the  digital  computer 
program  which  was  used  for  the  complete  Hybrid  computer 
simulation  and  controller  implementation.   The  program  includes 
all  parts  described  in  the  earlier  sections  and  appendixes. 
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